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Q^. We study single-spin asymmetries in semi-inclusive deep inelastic scattering with transversely 
(D , 

Z/J \ polarized target. Based on the QCD factorization approach, we consider Sivers and Collins contri- 
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butions to the asymmetries. We fit simple parameterizations for the Sivers and Collins functions 



■ to the recent HERMES data, and compare to results from COMPASS. Using the fitted parame- 
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. terizations for the Sivers functions, we predict the single transverse spin asymmetries for various 

CN ' processes in pp collisions at RHIC, including the Drell-Yan process and angular correlations in di-jet 

■ and jet-plus-photon production. These asymmetries are found to be sizable at forward rapidities. 
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I. INTRODUCTION 



Single-transverse spin asymmetries (SSA) in hadronic processes have a long history, start- 
ing from the 1970s and 1980s when surprisingly large SSAs were observed in p^p —>■ nX ^ 
and pp — > A|X [2j at forward rapidities of the produced hadron. They have again attracted 
much interest in recent years from both experimental and theoretical sides P]. In particular, 
first measurements bv the STAR, PHENIX, and BRAHMS collaborations at RHIC have now 
become ava„aHe BIQ. wM.K expend .Ke SSA ob.™«o„. fro. ..e flxed.a,. e„.«. 
range to the collider regime. Again, large asymmetries were found in p^p —>■ ttX at forward 
rapidities of the produced pion. Meanwhile, experimental studies in Deep Inelastic Scatter- 
ing (DIS) by the HERMES collaboration at DESY, by SMC at CERN, and by CLAS at 
the Jefferson Laboratorvalso show remarkably large SSAs in semi-inclusive hadron produc- 
tion, 7*p| ttX 0, y, 9 ! Iiol 111 I ■ Data from COMPASS for scattering off deuterons have 
been published as well l^f, which show no large asymmetry. On the theoretical side, there 
are several approaches to understanding SSAs within Quantum Chromodynamics (QCD) 
P, Q, . Recent interest focuses on the role of partonic transverse momentum in creating 
the observed asymmetries. Transverse-momentum-dependent (TMD) parton distributions 
and fragmentation functions, and their relevance for semi-inclusive DIS (SIDIS), the Drell- 
Yan process, and single- inclusive hadron production at hadron colliders have been investi- 
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27l |. Compared to the normal integrated 



distributions, the TMD distributions provide much more information; for example, some of 
them contain information on orbital angular momenta of partons in the nucleon and have 



also been linked to spatial distributions of partons l2^ 

n 

The Sivers function is one of these interesting TMD parton distributions. It represents 
a distribution of unpolarized quarks in a transversely polarized nucleon, through a correlation 
between the quark's transverse momentum k± and the nucleon polarization vector Sj_. The 
existence of the Sivers function requires final/initial-state interactions, and an interference 
between different helicity Fock states of the nucleon. In the absence of interactions, the Sivers 
function would vanish by time-reversal invariance of QCD, hence it is often referred to as a 
"naively time-reversal-odd" distribution. As was shown in 2J, |25|, |26| , the interactions are 
represented in a natural way by the gauge link that is required for a gauge-invariant definition 
of a TMD parton distribution. Interference between different helicity Fock states implies 
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nonzero orbital angular momentum [2j, |28[. Both these properties motivate the study of 
this function. The Sivers function will contribute to the target SSA in semi-inclusive DIS, 
but also to SSAs in polarized pp scattering processes such as the Drell-Yan process and di-jet 
or jet-photon correlations. We will discuss all these asymmetries in this paper. 

The Collins function is another "naively time-reversal-odd" function. It is a transverse- 
momentum dependent fragmentation function and was introduced in [2^. It represents a 
correlation between the transverse spin of a fragmenting quark and the transverse momen- 
tum of the hadron relative to the "jet axis" in the fragmentation process. Like the Sivers 
function, it vanishes when integrated over all transverse momentum. Indications of a non- 
vanishing Collins effect have been found in semi- inclusive DIS j9|. Very recently results for 
measurements in e~^e~ annihilation to two hadrons have been reported, which give direct 
evidence for the Collins effect js^. 

The formulation and study of TMD functions is really useful only when they appear in 
processes for which QCD factorization at small transverse momentum can be established. 
The processes, therefore, also need to be characterized by a large momentum scale, and 
there has to be additionally a small measured transverse momentum. Rigorous theoretical 
analyses of such reactions started from Collins and Soper's seminal paper in which 
they proved factorization for di-hadron semi-inclusive processes in e~^e~ annihilation. Non- 
perturbative TMD fragmentation functions were defined and then further studied along 
with TMD parton distributions in The approach was extended to Drell-Yan dimuon 
production at hadron colliders 

More recently, these factorization theorems and the one for semi-inclusive DIS have been 
carefully (re-)examined in the context of the gauge- i nvariant definitions of the TMD par- 



ton distributions and fragmentation functions 
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32| . paying in particular attention to the 



"naively time-reversal-odd" functions. In summary, QCD factorization has been established 
for three classes of semi-inclusive processes: di-hadron production in e~^e~ annihilation, semi- 
inclusive DIS, and the Drell-Yan process. It still remains to be seen whether factorization 
holds for more complicated processes in hadronic scattering, such as for di-jet (di-hadron) 
azimuthal angular correlations js^l^]. These processes, too, are characterized by a large 
scale (the individual jet or hadron transverse momenta), and by an additional small trans- 
verse momentum related, for example, to the pair transverse momentum, or to the deviation 
of the two jets from being "back-to-back". Note that this is in contrast to single-inclusive 



3 



processes at hadron colliders like p-^p irX. The spin asymmetries for such reactions are 
power-suppressed ("higher twist"), and the theoretical description should be based on the 
methods developed in where factorization in terms of higher- twist correlation functions 
was established. 

In this paper, we will use the factorization approaches at small transverse momentum 
discussed above to study the single spin asymmetries in semi-inclusive DIS. We will focus on 
the Sivers and Collins asymmetries which presently are the most interesting ones. We realize 
that at the current stage it is difficult to apply the full factorization formulas developed 
in the literature in fitting the data and making predictions. Instead, we will make some 
approximations, neglecting higher order terms in the hard and soft factors. In this way we 
of course introduce some theoretical shortcoming, which we hope can be overcome in future 
studies. Our purpose is to look for the "main effects", that is, to provide a quantitative 
description of the spin effects now studied experimentally, and to draw our conclusions on the 
Collins and Sivers functions from these. Another goal of this paper is to use the information 
gathered from semi-inclusive DIS to make predictions for processes at RHIC, which is now 
taking data in transversely polarized pp scattering. With the fitted parameterizations for 
the Sivers functions that describe the HERMES data very well, we will calculate the Sivers 
asymmetries for processes at RHIC, including Drell-Yan dimuon production and di-jet and 
jet-plus-photon correlations. We will demonstrate that these asymmetries are expected to 
be large at RHIC and should therefore be closely investigated in the future polarized pp 
runs. This would then provide further tests of the physical picture behind the SSA, and of 
our theoretical understanding. 

To predict the SSAs at RHIC from the distributions fitted in DIS relies on the factorization 
for the pp processes which, as discussed above, is so far only established for the Drell-Yan 
reaction. It also relies on the universality of the TMD distributions. This issue has been 
addres^d detaU ,n QHHSQ. It was found, fo. e.a,„p.e. tKat tKe S.ve. fu.ct.ons 
for the Drell-Yan process will have an opposite sign compared to those for SIDIS, as a result 
of the behavior of the gauge-links in the functions under the time-reversal operation. We 
will use this additional sign in our prediction for the Drell-Yan SSA at RHIC based on our 

n 

Sivers function fitted to the DIS data. However, recent work has shown that the issue of 
universality appears to be much more complicated for the case of di-jet correlations, where 
the more involved color structure has profound consequences on the gauge links. As a result. 
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the Sivers functions for this reaction will differ from those in DIS by more than just a sign. 
This reservation notwithstanding, in order to obtain an order of magnitude estimate we will 
assume in this paper that the Sivers functions to be used for di-jet correlations have the 
same sizes as those for the Drell-Yan processes, and opposite signs with respect to the DIS 
Sivers functions. 

In all calculations of cross sections and asymmetries below, we will use the GRV LO 
parameterizations for the unpolarized quark distributions 35| , and the Kretzer set of unpo- 
larized quark fragmentation functions 3]- These will also serve as starting points for our 
parameterizations of the Sivers and Collins functions. 

The rest of the paper is organized as follows. In Sec. II, we will review the basic formulas 
for the SSAs in SIDIS, and make model parameterizations for the Sivers and Collins func- 
tions. We then fit our parameterizations to the HERMES data. We will also compare our 
fit with the recent COMPASS data on the Sivers and Collins asymmetries. In Sec. Ill, we 
will calculate the Sivers asymmetries for the Drell-Yan process and for di-jet and jet-photon 
correlations at RHIC, using the fitted parameterizations from Sec. II. We summarize in 
Sec. IV. 



II. SSA IN SEMI-INCLUSIVE DEEP INELASTIC SCATTERING 

In this section, we will study the SSA in the SIDIS processes ep — > ehX and fid — > fihX, 
where h represents a hadron observed in the final state. We will compare the theoretical 
calculations of the asymmetries with the HERMES measurements. We will use some simple 
parameterizations for the Sivers functions and the Collins fragmentation functions, and fit 
these to the experimental data. A comparison of our fit with the COMPASS measurements 
will also be presented. Similar phenomenological studies of these asymmetries have also 
been performed in for the Sivers case and in ,39| for the Collins asymmetry, using 

the earlier HERMES data. 

We will start by briefly recalling the factorization formulas for the SIDIS process. For 
details, we refer the reader to Ref . j3l| . As discussed in the introduction, we will make some 
simplifying approximations, in order to sharpen the constraints on the Sivers and Collins 
functions. 
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A. Theoretical Formalism and Approximations 



The differential cross section for SIDIS, including the unpolarized part and the Sivers 
and Collins asymmetry contributions, may be written in the following form: 

da ^vra^jj^ 
dxBdydzhd'^Ph± 

-{l-y)xB\S^\sm{cPh + cPs)F^fH , (1) 



[l-y + yy2)xB [Fuu - sm{^h - 05)15^x1^^^") 



where (ph (05) is the angle between the lepton plane and the 7*-hadron-plane (and the 
transverse target spin), y is the fraction of the incident lepton energy carried by the photon, 
and Ph± is the (measured) transverse momentum of the hadron. In order to compare with 
the experimental data, in the above formula and the following calculations, the azimuthal 
angles (05 and (ph) are defined in the so-call virtual photon frame where the virtual photon 
is moving in the z direction. These definitions are different from those in where a 
hadron frame has been chosen to define these an gles . This difference has led to the different 
signs in the above formula, compared to that in j3l|. The structure functions Fm; and Fut 
will depend in general on Phj_, and on the invariant mass of the virtual photon, the 
Bjorken variable xb, and on the fraction of the photon longitudinal momentum carried 
by the hadron observed in the final state. According to the factorization formula of , the 
structure functions can be factorized into the TMD parton distributions and fragmentation 
functions, and soft and hard factors. For example, for the unpolarized structure function, 
we will have j^| 



FuuixB,Zh,Q'^,Ph±) = '^l j d^k^d^p^d^X 



q=u,d,s, 

xq {xB,k±,fi'^,XB(,p) q (^Zh,p±, fi^,C/zh, S{X±,fj,^,p) 
xH {Q^Ip\p) 6^^\zhk^+p^ + Ax - Phi.) . (2) 

This form is valid at low transverse momentum Phi_ ^ Q and is accurate at the leading 
power of -P^_|_/Q^. As seen from the 5-function expressing transverse-momentum conserva- 
tion, the observed hadron's transverse momentum is generated by three contributions: the 
transverse momentum k_[_ of partons in the nucleon, (described by the TMD distribution g), 
the transverse momentum px acquired in the fragmentation process (as expressed by the 
TMD fragmentation function g), and the combined transverse momenta A_l of (large-angle) 



soft-gluon radiation, embodied in the soft factor S. Each of these transverse momenta is in- 
tegrated in Eq. (0), but leaves its imprint in the distribution in Ph± of the observed hadron. 
In contrast to the TMD functions in (j2)), H is a hard factor that depends solely on the 
large scale Q. Furthermore, /i ~ Q is a renormalization scale, p a gluon rapidity cut-off 
parameter, and ( is defined as ( = {2P ■v)'^/v^, with P the target hadron momentum, taken 
in the "plus" -light-cone direction, P = {P^,P^ = 0, P± = 0±), and v a time-like vector 
conjugate to P, i.e., with only a "minus" -light-cone component. For details, see 3l|, where 
also the related definition of ( is given. Similar factorization formulas as (j21) can be written 
down for the single-transversely polarized structure functions F^*^'^^'* and F^'^'*"^. 

For simplicity, in the following numerical calculations, we will use the leading order 
expressions for the hard scattering and the soft factors, for which we have S = H = 1. 
At this order, we may also neglect the ( dependences in the parton distributions and 
fragmentation functions. All this brings us to the parton model picture for semi-inclusive 

n 

DIS [23]. However, we stress that higher order effects can be systematically and consistently 
studied only within the complete factorization framework. With the above approximations, 
the structure functions can be simplified to Isij ] 

Fuu = J d'^k±(fp_i_q{xB,k±)q{zh,p±)S^'^\zhk±+p±-Ph±) , 
F^r" = / d'k^d^P^^^^^qr {xb, k^) q {zh,p±) S^^\zhki_ + p^ - Ph±) , 

F^f'^^ = J d^k^d^pj_^^j^6qT (xB, k^) 6q 6^^\zHkj_ + p^~ Phi) , (3) 

where a sum over all quark and antiquark flavors, weighted with the squared quark electric 
charge, is implicitly understood from now on. Phi_ denotes a unit vector in direction of 
PhL- In the above equations, q and q represent the unpolarized quark distribution and 
fragmentation functions, respectively, the Sivers functions, 5q the Collins functions, and 
5qT the transversity distribution functions. The definitions of the above distributions and 
fragmentation functions are consistent with the so-called "Trento conventions" ko|, while 
opposite in sign with respect to that used in for the Sivers function. To optimize 
statistics, the experimental measurements of the asymmetries are normally presented after 
integrating over the modulus of the hadron's transverse momentum Ph^- After integration. 
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the cross section can be written as 



da da juj . , . daf}^^'^ . , , , n daf?^^^^ , . 

- sin(0/, - 05) -^-tS sin(0/, + 0s) ^ ^ , (4) 



dxBdydzhd(j)h dxsdydzh dxBdydzh dxsdydzh 

where the various differential cross sections will depend on x^, Zh and y. The dependence 
of the cross sections on the azimuthal angles results in the azimuthal asymmetries measured 
in experiment. The unpolarized cross section is given by 

dauu 47ra2^s / 



1 - 1/ + ^ XBq{xB)q{zh) , (5) 



dxBdydzh Q"^ \ 2 

where q^xb) and q{zh) are the integrated parton distribution and fragmentation functions. 
Here we assume that we can obtain the integrated parton distribution by integrating over 
the transverse momentum in the corresponding TMD parton distribution. This assumption 
will of course need to be modified if higher-order corrections are considered [SI]- Similarly, 
we can calculate the polarized cross sections. In these calculations, we further assume that 
the final hadron's transverse momentum is entirely related to the transverse-momentum 
dependence in the Sivers and Collins functions. The transverse momentum contributed by 
the other factors in the factorized formula will give some smearing effects which may be 
viewed as "sub-dominant" . After this approximation, we can write down the polarized cross 
sections as 

J ^sivers /|_„,2 „ „,2 

d^^d^ - '^^'^^^^~^ + y^^^^^ (^B)g(./.), (6) 

J Collins A^n"^ q 

-^.IS.ll^d (7) 

where 6qT is the integrated transversity distribution function. q^^'^\xB) and 6q''^/'^\zh) are 
defined as 



5q'^"'\zH) = j dW-^qT{z,,p^) . (8) 

The above formulas (jH)-® will be used in the following calculations to study the experi- 
mental data for the asymmetries as functions of and Zh- Before doing so, we need to set 
up models for the Collins and Sivers functions. 

We would like to add one more comment before we proceed. In the derivation of Eq. ® , 
we have omitted the transverse momentum dependence of the fragmentation function, which 
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we referred to as "sub-dominant". This "sub-dominant" contribution could become impor- 
tant at small Zh where we cannot neglect the influence of the transverse momentum in the 
fragmentation process. For example, for a typical transverse momentum of the flnal state 
hadron of Ph± ^ 200 ~ 300 MeV, the quark transverse momentum could be as large as 
2 ~ 3 GeV at ~ 0.1, much bigger than the typical value of the intrinsic quark transverse 
momentum for the Sivers function, which is of order of a few hundred MeV. This means 
that at small Zh our approximation will break down, and the transverse momentum in the 
fragmentation function will be important. This effect will smear out the polarized cross sec- 
tion and suppress the asymmetry. Indeed, when assuming a Gaussian transverse momentum 
dependence for both the distribution and the fragmentation functions, an additional factor 



of Zh appears, suppressing the polarized cross section at small Zh |38|. As a consequence, we 
should be cautious to apply Eq. (jH)) at small Zh- On the other hand, the above drawback 
does not apply to the case of the polarized cross section for the Collins contribution, Eq. (|7j), 
where we omit the transverse momentum dependence in the distribution. This is because 
there is no kinematic enhancement associated with the intrinsic transverse momentum in 
the parton distribution, compared to the fragmentation case. 

B. Model for the Sivers functions 

There exist by now quite a few model calculations for the quark Sivers functions in the 
nucleon The results of these vary rather widely. Here, we will instead adopt simple 

parameterizations for the Sivers functions and flt these to the HERMES data. We choose a 
form that has only a single free parameter for each flavor; the present data probably do not 
yet warrant a more complex form. Our parameterization is as follows: 



^^^(l - x) , (9) 



u{x) 

where in both equations u{x) is the unpolarized u-quark distribution. We assume that only 
the quark Sivers functions are non-zero, and that the antiquark ones vanish. This assumption 
will of course likely need to be modifled at small x. In the above parameterization, the factor 
X on the right-hand-side represents the valence nature of the Sivers function, whereas the 
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factor {1 — x) denotes an expected suppressed behavior of the function at large x ^ 



C. Model for the Collins functions 

For the Colhns asymmetry, we have two sets of unknown functions: the transversity dis- 
tributions and the Colhns fragmentation functions. From the present experimental data on 
the Collins asymmetries we cannot obtain constraints on both of them simultaneously ^. 
Here we adopt a parameterization for the transversity function [3] that is based on satu- 
ration of the Soffer inequality 4^. We note that this parameterization represents an upper 
bound for the transversity functions for the quarks. 

As for the Sivers functions, there have also been several model calculations for the Collins 
functions 4^, showing rather wide variations. Again, we will just use a simple parameter- 
ization for the Zh dependence of Sq^^^'^\ The flavor dependence of the Collins functions is 
important since one would like to describe the asymmetries for different hadron species. 
From the theory side, one could get constraints for the flavor dependence based on momen- 
tum conservation in the fragmentation process: the Schafer-Teryaev sum rules These 
sum rules state that the odd-moment (for the intrinsic transverse momentum) of the Collins 
function vanishes when the function is summed over all hadron states. Because the sum 
rules only involve the integrals over and p± (with a weight pj_) of the Collins functions, 
one cannot obtain from them more detailed constraints, e.g., for the Zh or pj_ dependences. 
In the following we will motivate a simple conjecture for the Collins functions, based on 
quark-hadron duality in the fragmentation process. This will provide us with additional 
constraints. The main result is that any quark Collins fragmentation function is very small 
when summed over all final hadrons. For example, the u quark Collins functions to all 
hadron final states will satisfy 

$^5M"(^^,P±) = 0(m„) . (10) 



^ This power suppression in (1 — x) could actually be as strong as (1 — x)"^ [33 ■ In our fit, most data are in 
the intermediate range of x, and the a; — s- 1 limit is not really reached. In any case, the power of (1 — x) 
will be modified by logarithms. 

^ Note, however, that independent information on the Collins functions is now coming from measurements 
of hadron-pair production by the BELLE collaboration • It is hoped that combination of these results 
with those from lepton scattering would eventually give information on transversity. 
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The above equation is motivated as follows. The Collins function is defined as 

X (o|/:(-oo, 0)^(0) |p,x) (P,x|^(e+, 6)/:t(^+, _oo) |o) 

where fc^ = / Zh and fc_L = —'p±/zh- Here, the final state hadron has been taken to have a 
large light-cone "minus" momentum component [we remind the reader that Ph denotes 
the momentum of the observed hadron, while is the (integrated) transverse momentum 
that the hadron acquires in the fragmentation process relative to the fragementing parton; 
see Eq. ©j. Finally, C is the gauge link along the light-cone direction conjugate to Ph- 
Here, Now, if we sum over all hadrons in the above Collins fragmentation functions, the 
intermediate hadronic states can be replaced by a quark or a quark plus gluons (or quark- 
antiquark pairs) using quark-hadron duality arguments. We then get the following equation: 

y'^^A^qh^^ ) ^ !L / ^/Le-('=-«+-^^-'^x)yTr{757-y (12) 

h X 

X (0|/:(-oo, O)^(O)lPgX) {P,X\^{t. h)C\t. b, -oo) |0) 

whose validity rests on the argument of quark-hadron duality for the fragmentation process. 
Duality-breaking effects will somewhat modify the above equation. The right hand side 
of Eq. ()12|) may be viewed as a quark Collins fragmentation function into a quark (or 
antiquark/gluon) state. The helicity-fiip required for a non- vanishing Collins function is 
then possible because of a finite quark mass. Thus, we approximately expect 

J2stizh,P±)=Oim,)^0 . (13) 

h 

If we further assume that the fragmentation functions for u and d quarks to strange mesons 
are suppressed relative to those into pions, we can have even stronger constraints for the 
pion Collins functions: 

Sq''^{zh,P±) + Sq^-{zh,Pi.) + Sq^\zh,p±) ^ , (14) 

where q represents any flavor of u,d quarks and their antiquarks. Further simplification of 
the above equation can be derived by considering isospin and charge symmetry relations 
between the different fragmentation functions. For example, we will have the following 
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6u^^ 


= 5(t' 


= 5(r^ 


= 




= 5u^~ 


= 6u^^ 


= 5(F 


6u-' 


= 


= 5d-' - 





relations, 



^unfav. 

^Mav.+S^nfav.] ■ (15) 

Here SqJ^^ and Sq^^j^^ represent the "favored" (in the sense that the leading Fock state of 
the hadron contains the parent quark flavor) and "unfavored" (where it does not contain 
it) fragmentation functions, respectively. Substituting the above relations into Eq. (fT^ . we 
will obtain 

S<i%.. + S^lnfav. ^ , (16) 

which means that the unfavored Collins function is approximately equal to the favored 
one with opposite sign. This result is counter to the usual notion in the literature that 
the favored fragmentation functions should be much larger than the unfavored ones. It of 
course crucially depends on the validity of the approximations made in the above derivation 
and will be subject to some corrections. We note, however, that this observation has also 
support from a string model description for the Collins fragmentation function As a 
test, we will treat in the following the favored and unfavored Collins functions free from the 
constraint (jlfij) . and fit them to the data to see if they naturally satisfy the above relation 
or not. To parameterize the Collins functions, we use the following two sets of functional 
forms. 

Set I : Sf/^J'\z) = Cfz{l - z)u-\z) , 5d2{z) = C^z{l - z)u-\z) , 
Set II : )(z) = Cfz{l - z)r\z) , 5gf/i(^) = Cuz{l - z)d-\z) . (17) 

The z factor in these parameterizations represents the vanishing of the Collins function at 
small z, and the (1 — z) factor follows arguments made in [2^. The difference between these 
two sets is that for Set I we parameterize both favored and unfavored Collins functions in 
terms of the favored unpolarized quark fragmentation function, while for Set II we parame- 
terize the unfavored Collins function using the unfavored unpolarized quark fragmentation 
function. Set I is inspired by the constraint of Eq. (jTHI); note that this ansatz is expected 
to violate the positivity constraints at very large Zh. On the other hand. Set II respects the 
positivity constraints, provided \Cf^y\ < 4. In the following, we will fit the HERMES data 
with these two sets of parameterizations. 
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D. Comparison with SIDIS Data 



We will now calculate the Collins and Sivers asymmetries using the functions specified 
above, and fit the free parameters to the new HERMES data on the Collins and Sivers 
asymmetries Here we will use = 2.41 GeV^, which is the average for the HERMES 
kinematics. We choose fi = Q in the unpolarized parton distribution and fragmentation func- 
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tions. We will then compare our fit results to the recent measurements by COMPASS 
For the fitting, we use the CERNLIB MINUIT routine H- 

As a function of a;^, the Sivers asymmetry can be calculated from the following formula: 

/ dzhdy- 



J II a dXBCLydZfi 

where the minus sign results from the minus sign in the polarized differential cross section in 
Eq. (jH). Since the y integral is the same for the numerator and denominator, it cancels out. 
Moreover, the integral over Zh can be factored out as a consequence of the approximations 
that led to Eqs. (0) and ® , and the Sivers asymmetry will be proportional just to the ratio 
of the Sivers functions over the unpolarized quark distributions, summed appropriately over 
fiavors. 



The Sivers asymmetry as a function of Zh can be calculated similarly. In Fig. d we s 
the results of our fit of the Sivers asymmetries for 7r+ and tt^ to the HERMES data 
For the two free parameters the fit gives 



low 
l3- 



5„ = -0.81 ± 0.07, 5rf = 1.86 ±0.28 , (19) 

with x^/d.o.f ~ 1.2. The band in each plot of Fig. 1 corresponds to a l-a error in the 
determined parameters. One can see that the u quark Sivers function appears somewhat 
better constrained by the data than the d quark one. This is readily understood from 
the fact that u quarks in DIS enter with the charge factor 4/9, so that in scattering off a 
proton target the u quark distribution is particularly selected. Another feature is that the d 
Sivers function comes out larger (by about a factor of two) than the u quark one, and with 
opposite sign. This behavior is quite different from model calculations [41]. The result is 
due to the fact that the HERMES Sivers asymmetry for vr" is much smaller than that for 
7r+. Theoretically, however, tt^ production should also have a significant contribution from 
u quarks, because one finds ~ O.Gm'^^ for the fragmentation functions when integrated 
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FIG. 1: Sivers SSA fit to the HERMES data flOi l: see text. The bands correspond to the 1-a error of 
the fitted parameters. Note that the data have not yet been corrected for acceptance and smearing. 

over the experimentally relevant region 0.2 < z < 0.7. To obtain a much smaller asymmetry 
for TT^ than for vr"^, there then have to be fairly strong cancellations between the u and d 
quark Sivers functions. We note that the signs we find for our Sivers functions are consistent 
with expectations in [2^, where they were qualitatively related to the opposites of the quark 
contributions to the proton anomalous magnetic moment. 

Figure 121 shows predictions for the vr*^ Sivers asymmetries as functions of xb and z^, based 
on our fits for the Sivers functions. We note that our prediction for the vr" asymmetry 
is nearly independent of Zh. This is because the u and d quark fragmentation functions 
for 7r° are the same, and because in our approximation the distribution and fragmentation 
functions are decoupled. However, we have to keep in mind that this decoupling might break 
down at small Zh, as we discussed before. This could be tested by future HERMES data. 



In m 



38l |. earlier HERMES results |9tl for the Sivers asymmetries were fitted. The meth- 



ods somewhat differed from ours. In 



a particular transverse momentum dependence is 



assumed for the Sivers functions and the unpolarized quark distribution and fragmentation 
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FIG. 2: Predicted Sivers SSA asymmetries for production at HERMES. 



functions. With more free parameters the experimental data are fitted equally well, and the 
u and d valence Sivers functions are obtained from the fit. In the asymmetries weighted 
with the transverse momentum of the hadron were used for the fit. Both fits find a large d 
quark Sivers function with opposite sign relative to the u-quark one. We have also checked 
that our fit results for the Sivers functions are consistent with these fits within the current 
large uncertainties, where we notice that the Sivers function in js^ has an opposite sign 
compared to ours and to the "Trento conventions" 1401 1. 

n 

The COMPASS collaboration also has measured the Sivers asymmetry J12l], separately 
for positively and negatively charged hadrons, produced off a deuteron target. To simplify 
the comparison with their data, we assume that the leading hadrons are mostly pions. We 
calculate the Sivers asymmetries for 7i~^ and vr" in the kinematic region of the COMPASS 
experiment, using the above fitted Sivers functions for u and d quarks, and compare to 
their data for leading positive and negative hadrons, respectively. We show this comparison 
in Fig. El One can see that our calculations based on fits to the HERMES data are also 
consistent with the COMPASS data, within error bars. We note that for the kinematical 
region of the COMPASS experiment, our predicted Sivers asymmetries for a deuteron target 
are very small, except in the large-x valence region. The smallness of the Sivers asymmetry 
is again related to cancellations between u and d contributions, which for deuterons enter 
in a different combination than for a proton target. It will be very interesting to check 
these predictions with future COMPASS data for a proton target. Thanks to the higher Q^, 
such data would also help in confirming the leading-twist nature of the Sivers and Collins 
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FIG. 3: Sivers asymmetries compared to the COMPASS data 



asymmetries. 

We next turn to the Collins asymmetry. Here we follow a similar procedure as we did for 
the Sivers case above. As we mentioned earlier, the situation is more complicated because of 
the fact the nucleon transversity densities are currently not known, and we need to resort to 
a model or ansatz for the latter. As described above, we will use the parameterizations for 
the quark transversity distributions of ■42], which represent upper bounds for the densities. 
We will fit to the HERMES data using the two sets of simple parameterizations for favored 
and unfavored Collins functions given in Eq. p7|) . 

The asymmetry as a function of is calculated from the formula 



q=u,d 



(20) 



q=u,d,u,d 

where again the minus sign comes from the sign in the polarized differential cross section, 
Eq. (Q. and q^ represent the fragmentation functions integrated over the accessed 

region in Zh- Kinematic cuts impose a correlation between xb and y, and the integral over 
y will depend on xb- In the experimental analysis, the data for the Collins asymmetries 
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FIG. 4: Same as Fig.n\for the Collins asymmetries, using the Set I parameterization of the Collins 
functions. The data ll (\J and the theory curves are for the so-called lepton beam asymmetries. Note 
that the data have not yet been corrected for acceptance and smearing. 

are presented in two different ways. One is to give results in terms of the virtual-photon 
asymmetry, factoring out the term {1 — y)/{l — y + 2/^/2). The other way is to give the 
directly measured lepton-beam asymmetry. In our calculations, we follow the latter way. 
We neglect the contribution of longitudinal photons to the unpolarized cross section, which 
HERMES has considered in the analysis of the virtual-photon asymmetries In view of 
the overall uncertainties, this is a minor effect, as we have checked by comparing also to the 
virtual-photon asymmetries. From the fit to the lepton-beam asymmetry data, we get the 
two fit parameters as follows: 



Set I : Cf = -0.29 ± 0.04, C„ = 0.33 ± 0.04 
Set II : Cf = -0.29 ± 0.02, = 0.56 ± 0.07 



(21) 
(22) 



with x^/d.o.f. ~ 0.8(0.7) for the Set I and Set II parameterizations, respectively. The fit 
results are shown in Figs. |3] and compared to the HERMES data. Both fits are of the 
same quality. 
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FIG. 5: Same as Fig. ^ for the Set II parameterizations for the Collins functions. 



In Figs, ini and [71 we plot the fitted favored and unfavored Collins functions (times z) for 
Sets I and II respectively. Note that we multiply the favored ones by (—1) to compare their 
magnitudes. For comparison, we also show the corresponding unpolarized quark fragmen- 
tation functions 36]. It is evident that the two sets of Collins functions indeed both satisfy 
the positivity constraints. The equal quality of the fits obtained for sets I and II implies 
that the current experimental data neither necessarily support the constraints we derived in 
Eq. (fT^ . nor do they rule them out. However, from both fits we indeed find that in a quite 
large range of Zh the unfavored Collins function has the same size as that of the favored 
one with opposite sign. A similar conclusion was obtained from a fit to this asymmetry 
using the transversity functions calculated in the Chiral Quark Model [s^. We hope that 
higher-statistics data will become available in the near future that will test the relations. 

Figures IHl and El show predictions for the 7r° Collins asymmetries as functions of and Zh, 
based on our fits for the Collins functions. From these plots, we find that the asymmetries 
are very small for both sets of the Collins functions, because of strong cancellations between 
the contributions from favored and unfavored Collins functions. We note that preliminary 
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FIG. 6: Set I { — l)x favored and unfavored Collins fragmentation functions. Also shown are the 
unpolarized quark fragmentation functions from Kretzer's parameterizations. 
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FIG. 7: Same as Fig. \^for the Set II Collins fragmentation functions. 

data from HERMES indeed indicate that the 7r° asymmetry is consistent with zero 
Higher-statistics data on this will be highly interesting for testing the above conclusions. 

As before, we also compare our fit to the COMPASS measurements [12] for the Collins 
asymmetries, where only the virtual-photon asymmetries are presented. We also note that 
the convention for the Collins asymmetry used by the COMPASS Collaboration is different 
from that used by HERMES. In our calculations, we have made the relevant modifications 
in order to compare with the COMPASS data. We show these comparisons in Figs. ^1 
and ^2 for our Set I and II Collins functions, respectively. As for the Sivers case above, 
there is good consistency. The overall asymmetries are again small because of the deuteron 
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FIG. 8: Predicted Collins SSA asymmetries for vr^ production at HERMES with Set I. 
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FIG. 9: Predicted Collins SSA asymmetries for tt^ production at HERMES with Set U. 

target used, and because the assumed transversity sea quark distributions are small. Future 
high-statistics COMPASS data for a proton target would be highly interesting. 

III. SINGLE-TRANSVERSE SPIN ASYMMETRIES AT HADRON COLLIDERS 



An important issue in the study of hard-scattering processes is the universality of the 
nonperturbative objects: the parton distributions and fragmentation functions. In the case 
of single-spin asymmetries, if universality holds, the Sivers functions obtained from, for 
example, SIDIS can be used to predict single-spin asymmetries in pp or pp scattering. The 
universality between different classes of processes is a complicated and interesting issue that 



has attracted much interest recently 



3^ 



49 



50|. As we described earlier, it was found 
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Qy for Set I q 



»/ our fitted Collins 



that, while strict universahty is violated already when going from SIDIS to the Drell-Yan 
process, the Sivers functions for the two processes differ only by a sign. It is therefore 
possible to use the fitted Sivers function from the last section and to predict the Sivers 
single-spin asymmetry for the Drell-Yan process at hadron colliders. More complicated 
processes in hadronic scattering, such as the SSA in di-jet angular correlations j^, are not 
yet completely understood at present, as far as factorization and universality are concerned. 
Progress has been made recently [3^; it appears that similarly to the Drell-Yan process 
universality is violated only by terms that are calculable from the color structure of the 
partonic scattering and hence may be taken into account in phenomenological analyses. 
Below, we will also give estimates for Sivers contributions to SSAs in di-jet correlations and 
in jet-plus-photon correlations at RHIC, using the Sivers functions of Section II, and the 
usual unpolarized hard-scattering functions. In the light of Ref. 3^, we expect that our 
estimates will likely need to be revised once these reactions will be completely understood 
in the context of factorization and universality, to take into account the appropriate factors 
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FIG. 11: Same as Fig. \1(A but for Set II of Collins fragmentation functions. 

embodying the nonuniversality of the Sivers functions. We shall briefly return to this point 
below. 

A. Drell-Yan Dimuon Production p^p —>■ fi^i^i^X 



In this subsection, we will calculate the Sivers single-spin asymmetry for the Drell-Yan 
process at RHIC, using the fit result of the last section [see Eqs. Q and (HD)]. As just 
discussed, one has 

q?^ = -q?^' . (23) 

After integrating out the lepton angles in the rest frame of the virtual photon, we obtain 
the following differential cross section for the Drell-Yan process: 

= ^ I'^'"'"" ^ ■ 

where M is the invariant mass of the lepton pair, q± the virtual photon's transverse mo- 
mentum, and y its rapidity. = 0^* — (ps is the difference between the azimuthal angles 
of the virtual photon and the transverse polarization vector in a frame where the polarized 
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hadron is moving in the z direction. At low transverse momentum, x\ and X2 are related to 
the mass and rapidity through x\ = M/ ^/s and X2 = M/ y/s where s is the hadronic 



center-of-mass energy squared. According to the factorization theorem j^, the hadronic 
tensors Wq and Wut can be factorized into the TMD parton distributions, and soft and 
hard-scattering factors. Again, neglecting the soft factor and using the Born expression for 
the hard part, we obtain simple expressions for the tensors: 

Wo = / d^hLd^k2±q{xi,kiA_)q{xi,k2A_)6^'^\ki^ + k2±- q±) , 

Wtu = J2jj ^'^iA±^^^^gT(xi,A;ix)g(xi,A;2x)5(^^ , (25) 

where qrixi, ki±) is now the Sivers function for the Drell-Yan process. A further approx- 
imation can be made by integrating out the transverse momentum \q±\, but keeping the 
dependence on azimuthal angle. The differential cross section can then be written as 

da Aa'^TT 
dM^dydcj) ~ 3sM^ 

where 

Wo = J2 y 5(^1)^(^2) , (27) 

q 

and 

W^f/T = 5^|gi^'/'^(xi)g(x2) . (28) 
<} 

In Fig. ^[ we plot the sin asymmetries as functions of the photon rapidity y and the invari- 
ant mass M. From this plot, we see that the Sivers SSA asymmetry for the Drell-Yan process 
at RHIC is expected to be sizable for large rapidity, and should be measurable at RHIC if 
enough statistics can be accumulated in transverse-spin running. We note that the Sivers 
asymmetry in the Drell-Yan process is also a particular focus for proposed measurements in 



Wo{xi,X2,M^) + sm(l)WTu{xi,X2,M^) , (26) 



polarized pp scattering at the planned GSI-FAIR facility |37l . l51 



52|. 



B. Correlations in p^p — > jeti(-Pi^) + jet2(-P2±) + X 

Other interesting observables at hadron colliders from which one can access the intrinsic 
transverse-momentum dependence of parton distributions are "back-to-back" correlations 
between two jets js^ Q. More specifically, we are interested in situations in which the 
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FIG. 12: Sivers asymmetries for the Drell-Yan process at RHIC, as functions of virtual-photon 
rapidity y and invariant mass M. 



sum of the two jet transverse momenta, q± = Pi± + P2± (or a component or projection 
thereof), is measured, while both Pij_ and individually are large. As for the Drell-Yan 
process discussed above, this is usually a small transverse momentum, much smaller than the 
large scales \Pi±\ ~ 1-^2x1 characterizing the overall process, and a special factorization may 
apply. Let us first, however, consider the cross section integrated over Here, coUinear 



factorization applies, and the di-jet cross section has the parton-model expression 



da 



da 



dyidy2dPl 



y^Xafa{Xa)Xbfb{Xb)—^{ab cd) 

^-^ dt 



ab 



(29) 



where da/dt is the differential cross section for the partonic process ab cd, with fa^^ 
the appropriate parton distribution functions. We have defined the transverse momentum 
P± = \Pi_l\ = I Ax 1 5 and yi and y2 denote the rapidities of the two jets. The kinematics are 
as follows: 



Xa 



+ e 



y2] 



Xb 



P± 



t 



1), 



u 



-Pi {e^'-y^ + 1) . 



(30) 



Here, s, i, and u are the usual partonic Mandelstam variables appearing in the partonic 
cross sections for the reactions ab — > cd, namely s = (Pa + Pb)^, t = {.Pa—PcYiU = {Pa—PdY, 
in obvious notation of the partonic momenta. The leading order contributions produce the 
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di-jet pair exactly balanced, that is back-to-back in the partonic center-of-mass frame. An 
imbalance in the transverse direction is generated by higher-order QCD corrections. At small 
but nonzero imbalance between the two jets, the dominant contributions will come from the 
intrinsic transverse momenta of the initial partons. As a model, we will generalize the above 
factorization formula to the case of small gl, in analogy with the SIDIS and Drell-Yan cases 
discussed earlier, taking into account the various contributions to the transverse momentum 
dependence coming from the parton distributions and soft factors: 
da 



dyidy2dP]^d'^qA 



^ j d'^ki±d'^k2±d'^X±XafaiXa,ki±)XbfbiXb,k2±) 
ab 

xSab^cdiX±)Hab^cdiPl)S^'\h± + h± + Ax - gl) , (31) 



where Sabred is a soft factor for the process ab — > cd, while Hab^cd is the hard part of the 
reaction, related to lowest order to da/dt. We emphasize the overly simplistic character of 
Eq. (jHT|) as it stands. The detailed factorized form (if it exists) will likely be different; in 
particular, one expects an interplay of the color structures of the soft factors and the hard 
parts, as found in resummation studies for jet cross sections Q]. 

In a similar fashion, we write the Sivers-type contribution daTu to the single-polarized 
cross section, 

J J = dauu + ez ■ (s± x q±) daru , (32) 

dyidy2dPld^q± V / 

where is the unit vector in the z-axis direction. In a factorized form, we will get: 
daru = y d'^h±d'^k2±d'^X±^^^j^XaqTaixa, ki±)xbfb{xb, k2±) 

ab 

Sab^cdiX±)Hab^cdiPl)5^'\h± + k2± + X± ' q±) . (33) 



X 



We can further simplify the polarized cross section by evaluating the expression 

I5i 



|Pj_| (sin 01 + sin 02) 



1^1 



\(1±\ 

^ \S^\ (^Sgn(7r - 9) cos0i + sm^^^^ , (34) 

where 0i and 02 are the azimuthal angles of the two jets relative to the polarization vector 
S±, and ^ = 02 — 01 the angle between the two jet transverse momenta. All these azimuthal 

25 



angles are defined in a frame tliat tlie polarized proton is moving the z direction. In the above 
derivation, we have used the approximations |-Pi±| ~ \P2\\ ~ and |g_L| ~ |Pi||sin6'|, 
which are valid at small g_L [Q is close to tt). From the above result, we can see that there 
are two terms contributing to the spin asymmetry: one is with cos 01 and the other with 
sin 01. The first term has a Sign function associated, which gives a positive contribution 
when Q is smaller than vr and a negative one otherwise. 

In the above formulas, we did not include any gluon Sivers function contributions. The 
gluon Sivers function could dominate the asymmetry at central rapidities js^. Another 
important issue is the relevant Sudakov suppressions for the asymmetries, which was found 
to be sizable for the di-jet correlation in the RHIC energy range Issl. In the following 
numerical studies, as an order of magnitude estimate for these asymmetries, we will neglect 
these effects, which however should be taken into account in future more detailed studies. 

Following the same procedure that we used for the Sivers asymmetries in SIDIS and 
Drell-Yan processes, we can further simplify the polarized cross section by integrating out the 
transverse momentum, but keeping the azimuthal angle dependence explicit. The differential 
cross section then can be written as 

27rc?cr ^ dayy ^ mz^^ ^^^^'^ I oin0, ^'^^^'^ (35) 
dyidy2dP]_d(f)i dyidy2dP]^ dyidy2dP]_ dyidy2dPl 

where the unpolarized cross section has been given in Eq. ()29|). and the polarized ones read 

da!^^ (u-?.) , s „ , .da 



^XaqTl^\xa)xbfb{xb)—^{ab cd) 



dy,dy2dPl ^ .-^a,-^...,-^., 



dy^dy2dPl \P^\ ^---lav^^,-^.^^^, ^£ 

ab 



with the distribution q^^ defined as 



4\x) = J d'k^^^qr{x,k^) . (37) 

We note that the second term in the polarized cross section is power-suppressed by M/P±. 
This suppression is due to the fact that we have integrated over all intrinsic transverse 
momentum. Clearly, this term is beyond the approximations we have made, and we cannot 
reliably predict it since there will be other sources of power-suppressed contributions, for 
example generated within the Qiu-Sterman mechanism Q|. Since it is anyway expected to 
be small, we will discard it in the following. Employing the same set of Sivers functions 
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that we used for our predictions for the Drell-Yan process above, we then find the results for 
An = da^lj / dauu shown in Fig. UHl We show the asymmetries as functions of the rapidity 
of jet "1", and of the jet transverse momentum P^. One can see that the SSA for the di-jet 
correlation can become very large, in particular in the forward rapidity region. Asymmetries 
of this size should be relatively easily measurable in the future. We also note that the 
asymmetry has opposite sign compared to that for Drell-Yan dimuon production discussed 
earlier. The reason for this is that w-quark contributions dominate in Drell-Yan, thanks to 
their large electromagnetic charge, whereas for di-jets, (i-quark contributions are not charge- 
suppressed and in fact dominate, keeping in mind that the analysis of the HERMES data 
appears to favor a large d quark Sivers function. The opposite signs of the Sivers up and 
down quark distributions we found in Eq. (jl9|) then explains the opposite signs of the spin 
asymmetries for Drell-Yan and di-jets. We note that if the two jets are within the central 
rapidity region, our prediction for the asymmetry is much smaller. As mentioned above, the 
gluon Sivers function could dominate the asymmetry in this region j^^. 

We stress again that if factorization can be shown for the Sivers SSA in di-jet production, 
it is likely that the structure of the resulting expression may differ from the one we use. In 
particular, there will be calculable factors that represent the non-universality of the Sivers 
functions related to the process dependence of the gauge links, leading effectively to modified 
partonic hard-scattering functions 13], at variance with our use of the standard unpolarized 
ones. As a test, we have also used the modified partonic cross sections derived in [s^ . 
We find relatively small changes in the results we obtain. Unfortunately, however, this is 
not really representative: the cross sections given in 2^ are only for the quark- (anti) quark 
scattering channels, whereas the dominant contribution in our calculation mostly comes 
from qg scattering (with the gluon from the unpolarized proton). It remains to be seen 
to what extent eventually our predictions will change, once the process-dependence for the 
Sivers functions in di-jet correlations is completely understood. 

We finally note that it would also be interesting - in particular for measurements with the 
PHENIX detector - to study correlations between hadrons in opposite jets. Such di-hadron 
correlations could serve as surrogates for the di-jet correlations we have discussed above. In 
this case there will, however, also be contributions from the Collins mechanism. 
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FIG. 13: Sivers asymmetries for di-jet correlations at RHIC, as functions of rapidity yi and trans- 
verse momentum P±. 

C. Jet-Photon Correlations in p^p — > jet + 7 + X 

It is straightforward to extend the analysis of di-jet correlations discussed above to the 
case of jet-plus-photon correlations. We simply need to implement the cross sections for the 
appropriate Born- level partonic scatterings qq — > jg and qg — > jq in Eq. ()36|) . Although 
events with a photon suffer from smaller rates than two-jet events, they would offer additional 
information on the Sivers functions. It is also likely that proofs of factorization are more 
easily obtained here, since the reactions qq 'jg and qg 7g each have only a single color 
structure. Figure IT^ shows results for the single-spin asymmetry for jet-photon correlations 
for the same kinematics as for the di-jet case in Fig. ^1 Variables with the subscript "1" 
denote photon variables. Again, sizable asymmetries are seen, in particular at forward 
rapidities. The asymmetries are somewhat smaller than the ones we found for di-jets. This 
is a result of cancellations between our Sivers u and d functions, due to the larger weighting 
factor 4/9 that the w-quark contributions now have for the prompt-photon case. 

IV. CONCLUSIONS 

In this paper we have studied single-transverse spin asymmetries in semi-inclusive deep 
inelastic scattering and at hadron colliders. We have analyzed the Sivers and Collins con- 
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FIG. 14: Sivers asymmetries for jet-photon correlations at RHIC, as functions of photon rapidity 
yi and transverse momentum P±. 

tributions to the spin asymmetry in SIDIS, and fitted simple parameterizations of tlie cor- 
responding functions to recent data from HERMES. These fits work well and also turn out 
to be consistent with COMPASS measurements of the asymmetries in DIS off a deuteron 
target. For the Sivers functions, we found dominance of the down-quark distribution over 
the up-quark one. The Sivers-ci density in SIDIS turns out to be positive, while the u-quark 
distribution comes out negative. Concerning the Collins functions, we have found that cur- 
rent data do not yet pin down the relative size of "favored" and "unfavored" functions, 
which is also due to the fact that the transversity densities are not yet known. We have also 
given theoretical arguments that the "favored" and "unfavored" Collins functions could be 
of similar size, and of opposite sign. 

We have then investigated Sivers-type single-spin asymmetries at hadron colliders, focus- 
ing on the Drell-Yan process and on di-jet and jet-photon correlations, all in circumstances 
where there is a small measured transverse momentum, but the process is overall character- 
ized by a large scale. Using the Sivers functions obtained from the analysis of the HERMES 
data, we have made predictions for single-spin asymmetries for these processes. We find rela- 
tively large asymmetries, in particular at forward rapidities of the observed final state. Such 
asymmetries should be measurable with dedicated efforts at RHIC. Besides the additional 
valuable information they would give on the Sivers functions and therefore on the structure 
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of the nucleoli, they would also provide a test of our theoretical understanding of "naively 
time-reversal-odd" phenomena in QCD. The crucial issues in this are the factorization of the 
corresponding cross sections, and the universality of the Sivers functions, on both of which 
further theoretical work is required. 
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Note Added: Upon completion of this paper, we noticed the preprint Q] where also 
the Sivers functions were fitted to the new HERMES and COMPASS data, and predictions 
for SSAs in the Drell-Yan process were made. As far as we can see, our results are in 



qualitative agreement with those of 



keeping in mind that their sign convention for the 



Sivers function is opposite to ours and to that in the "Trento conventions". 
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